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Abstract 

In this paper we want to extend the analysis of the thermodynamics of the chmate system to 
the investigation of the role of its various scales of variability. This is achieved by performing 
coarse graining of the climate fields that allows for separating the contributions emerging at the 
various scales. We show that the coarser is the graining, the lower is the resulting estimate of 
the material entropy production. In other terms, all the spatial and temporal scales of variability 
of the thermodynamic fields of the climate system provide a positive contribution to the material 
entropy production. This may be interpreted also as that, at all scales, the temperature fields 
and the material heating fields have a negative correlation, while the opposite holds between the 
temperature fields and the radiative heating fields. Moreover, we obtain the latter correlations are 
stronger, which confirms that radiation acts as primary drive for the fiuctuations in the thermody- 
namic properties of the climatic fields, while the material fluxes dampen such fluctuations through 
dissipative processes. We also show how, using speciflc coarse-graining procedures, one can sepa- 
rate the various contributions to the material entropy production coming from the dissipation of 
kinetic energy, the vertical sensible and latent heat fluxes, and the large scale horizontal fluxes, 
without resorting to the full three-dimensional time dependent fields. We find that most of the 
entropy production is associated to irreversible exchanges occurring along the vertical direction, 
and that neglecting the horizontal and time variability of the fields has a relatively small impact 
on the estimate of the material entropy production. These findings deserve to be explored in a 
more general context, because they provide a way for understanding forced fluctuations and dissi- 
pative processes in fluids at various spatial and temporal scales. Moreover, the approach presented 
here seems promising specifically for the investigation of the atmosphere of exoplanets, because it 
allows for evaluating the error in the estimate of their thermodynamical properties due to the lack 
of high-resolution data. 
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I. INTRODUCTION 

Developing a comprehensive theory of chmate dynamics is one of the grand contempo- 
rary scientific challenges, both because of its obvious environmental, social and economical 
relevance, and because of the fundamental difficulties resulting from the fact that the cli- 
mate is an outstanding example of a complex system featuring non-trivial dynamics on a 
vast range of spatial and temporal scales. In recent years, the extraordinary developments 
of planetary sciences coming from the discovery of extra-solar planets and the ensuing need 
for understanding the properties of atmospheric circulations realized under physical and 
chemical conditions very different from those of the Earth and of the other solar planets 
makes even more urgent the need for providing new ideas and new methods for understand- 
ing the dynamics of fluid systems in rotating frames of reference undergoing spatially and 
temporally inhomogeneous forcings. 

Along the lines of the powerful theoretical construction due to Lorenz [24, 25] of energy 
cycle of the atmosphere, the climate can be seen as a non-equilibrium system, which gen- 
erates entropy through a variety of irreversible processes [8, 38-40], and transforms moist 
static energy into mechanical energy, as it features a positive spatio-temporal correlation 
between heating and temperature patterns, so that it can be represented schematically as a 
heat engine with a given efficiency [14, 27]. For a given value of the external and internal 
parameters, the climate system achieves a steady state by balancing the input and output 
of energy and entropy with the surrounding environment [41]. Interestingly, the large scale 
motions of the geophysical fluids are at the same time the result of the mechanical work 
(continuously dissipated through the turbulent cascade) produced by the climatic engine, 
and contribute to reducing the temperature gradients which make the energy conversion 
possible [40, 45]. Obtaining a closure to this problem would be equivalent to developing a 
self-consistent theory of climate dynamics. 

Recently, a line has been drawn connecting the climate engine's efficiency, the intensity 
of the Lorenz energy cycle, the entropy production of the climate system. In particular, 
the efficiency of the equivalent thermal machine sets also the proportionality between the 
internal entropy fluctuation of the system and the lower bound to entropy production by 
the fluid compatible with the 2nd law of thermodynamics. Such a bound is basically given 
by the entropy produced by the dissipation of the mechanical energy, whereas the excess of 



entropy production is due to the transport of heat down the gradient of the temperature 
field. Therefore, it has been possible to introduce a index of irreversibility of the system 
able to characterize the relative importance of these two processes [27]. 

The tools of phenomenological non-equilibrium thermodynamics [2, 19] seem well suited 
for providing powerful tools for addressing the long-standing problem of characterizing the 
structural properties of the climate system [43]. The goal is to define a unifying framework 
for understanding the climate's variability, its feedbacks, and its large-scale processes, in- 
cluding the atmosphere-ocean coupling, the hydrological cycle, as well as understanding the 
mechanisms involved in climate phase transitions observed at the so-called tipping points, 
i.e. conditions under which catastrophes may occur for small variations in the boundary 
conditions or in the internal parameters of the system [22, 33]. While the thermodynamic 
interpretation of the baroclinic disturbances, which provide the dominant contributions to 
the low-to-high latitudes heat transport, lies at the core of dynamical meteorology [13], 
and thermodynamics provides indeed the best framework for studying strong meteorologi- 
cal features like hurricanes [6], recent results suggest that tipping points can be effectively 
analyzed using the thermodynamic indicators developed in [27], with the efficiency and the 
entropy production providing the most interesting indicators [1, 30]. Moreover, recent stud- 
ies have underlined that it is instead possible to define generalized climate sensitivities able 
to describe quite accurately the responses of thermodynamic quantities to changes in CO2 
concentration, and to link such responses to variations in the intensity of processes, such as 
the hydrological cycle, which are much more commonly considered in the context of climate 
change literature [28]. 

Despite its relevance at theoretical level, traditionally, entropy production is not one of 
the ffist physical quantities climate modelers investigate when assessing the performance of 
a global climate model (GCM), when intercomparing it to other GCMs, or when evaluating 
the impact of changes of the chemical or orbital parameters of the system, such as when 
studying anthropogenic or natural climate change on various temporal scales. Nonetheless, 
recently it has been shown that, starting from the theoretical formulations proposed in 
[8, 27, 37] ,it is possible to provide a slightly simplified expression for the entropy production 
based only upon the radiative balances at surface and at the top of the atmosphere. Such 
expression allows for separating the contribution to entropy production coming from vertical 
and horizontal irreversible exchange processes and can be used on one side for testing GCMs, 



and, on the other side, for estimating or constraining the entropy production (and other 
thermodynamical quantities) for planetary objects for which only coarse resolution data 
are available and model-assisted reconstruction is mostly temptative [29]. One must keep 
in mind that from estimates of the entropy production it is possible to construct bounds 
on the energetics of the circulation, which is a property of paramount importance for the 
definition of the characteristic of a planetary atmosphere [29] . 

One should note that the attention towards entropy production in the climate system and 
in climate modeling has been revived when several authors started proposing it as target 
function to maximize when tuning free or empirical parameters of approximate numerical 
models [17, 21, 26] or for getting good first order approximations of the climate state without 
resorting to long integrations [12, 33-35] This is the weak or pragmatic version of the so- 
called maximum entropy production principle [18], which, in its strong form, proposes that 
any non-equilibrium systems adjust itself in order to maximize the entropy production; 
see [18, 32] for rather comprehensive review of this literature. The theoretical foundations 
proposed for such a principle [3, 4] have been severely criticized both at theoretical level 
[10] and in terms of its geophysical applications [9, 37], so that weaker formulations are 
now mostly preferred [5], while at the same time there are attempts at using the entropy 
production of the system as an approximate Lyapunov functional for the dynamical system 
under investigation [11]. 

In this paper, we wish to investigate the entropy production of a climate model for 
studying, instead of large scale balances, its fluctuations at different temporal and spatial 
scales. Climate is a multi-scale system where dynamics takes place on vast range of inter- 
acting scales. The definition of parametrizations for unresolved scales is a major challenge 
of climate modeling and the proposal of closure theories connecting small and large scale 
properties is a major part of any attempt at formulating approximate theories for climate 
dynamics. The issue of understanding the impact of small scales on large scales and vice- 
versa, and of performing properly the upscaling and downscaling of a model's output is 
of great relevance also for intercomparing the performances of various versions of a given 
numerical model, or of a set of numerical model simulating the same system, differing for 
the adopted spatial and temporal resolution, and for comparing model data to observations. 

Our goal is manifold. One one side, we want to introduce a way to evaluate how the 
different scales of motion contribute to the overall entropy production of the climate system. 



This investigation, therefore, complements the investigation of how much energy is contained 
in the various scales of motion and of the energy fluxes across these scales. In order to 
achieve this goal, we consider the entropy budget of the FAMOUS atmosphere-ocean general 
circulation model[15] in standard, present climate configuration, taking advantage of the fact 
that it is one of the very few climate models where the entropy production diagnostics has 
been implemented and throughly tested [36, 37]. Starting from the output fields at the 
highest possible resolution given by the model (temporal resolution of one time step, and 
same spatial resolution of the actual numerical model), we perform a coarse graining in space 
and in time to the dynamical and thermodynamical fields appearing in the terms describing 
the entropy production of the system, and we test how the estimate of the entropy production 
changes when different coarse graining are applied to the data. We can anticipate that we 
obtain that the coarser is the graining procedure, the lower is the estimate of the entropy 
production one obtains, as somehow intuitive. One must note that this is, in fact an obvious 
result when considering simple diffusive system, but not so obvious when fully nonlinear, 
multiphase systems are considered. We also obtain similar results when considering different 
degrees of longitudinal averaging of the fields, up to considering zonally averaged fields 
only. Our findings provide a way to assess how having low-resolution information about 
the dynamics of turbulent systems affects our ability to reconstruct its thermodynamical 
properties. Moreover, the procedure discussed in this paper allows to put on firmer ground 
the results proposed in [29] on the possibility of separating vertical and horizontal exchange 
processes as far as entropy production is concerned. Finally, we can study in detail the 
relationship between two apparently equivalent ways of computing the entropy production 
proposed in [8]. 

This paper is structured as follows. In section II we briefly recapitulate some definitions 
and equations relevant for setting the problem of computing the entropy production of the 
climate system we explain how to perform such a calculation in a climate model. In section 
III we explain what we mean precisely by coarse graining of the data and describe how it 
is actually implemented in the model's output. We also provide some conjectures what will 
be discussed in later in the paper. In section IV we present our results. We first describe 
the impact of performing coarse graining on time alone, thus exploring the range between 
time step data and long term averaged data, and then we extend our analysis to the space 
domain, showing how performing zonal, horizontal, and mass-weighted averaging over the 
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output data impacts the obtained estimate of the entropy production. In section V we 
present our conclusions and perspective for future works. In appendix A we present some 
theoretical arguments on a simple diffusive system for clarifying the meaning of the results 
obtained from the data analysis. 

II. CLIMATE ENTROPY BUDGET 

Following [2, 19], for any system it is possible to decompose the rate of change of its 
entropy dS/dt as dS/dt = deS/dt + diS/dt, where the first term is called the external and 
the second term is the internal contribution to the entropy budget. The external contribution 
corresponds to the entropy flux through the boundaries of the system whereas the internal 
entropy production is associated with the irreversible processes taking place in the system. 
The second law of thermodynamics imposes that the internal entropy production has to be 
nonnegative at all instants, so that diS/dt > 0. When a statistically steady state is achieved, 
the external and internal entropy production have to balance each other so that the total 
rate of entropy change is zero: we have 



dS/dt = ^ diS/dt = -deS/dt > 0, 

where the overline indicates averaging over a long time interval compared to the internal 
scales of the system. The previous expression means that a non-equilibrium system generates 
on the average a positive amount of entropy through irreversible processes, and such excess 
of entropy is expelled at the boundaries. Non-equilibrium is maintained if the system is not 
in contact with just one reservoir at a given temperature and/or chemical potential [7]. Of 
course, if the system is at equilibrium, the previous inequality becomes an equality, as in 
the long run the system reaches an homogeneous state of maximum entropy and entropy is 
generated no more. 

In the climate system two rather different set of processes contribute to the total entropy 
production [8, 40]. The first set of processes are responsible for the irreversible thermali- 
sation of the photons emitted near the Sun's corona at roughly 5700 K at the much lower 
temperatures, typical of the Earth's climate. This contributes for about 95% of the total 
average rate of entropy production for our planet, which is about 900 W m~'^ [8, 40]. The 
remaining contribution is due to the processes responsible for mixing and diffusion inside 



the fluid component of tlie Eartli system, and for the dissipation of kinetic energy due to 
viscosity. This constitutes the so-called material entropy production, and is considered to be 
the entropy related quantity of main interest as far as the properties of the climate system 
are concerned. See [37] for an extensive discussion of this issue a careful estimate of its value 
in two climate models, including the one used in this study. 

When separating the entropy budget for radiation and for the fluid part of the climate 
system, and taking long term averages, one can derive the following equation [8, 14]: 



V 



Qrad 
rp I "I" ^mat 



(1) 



where the integral is over the whole volume V of the climate system, g^ad is the radiative 
heating rate, Smat is the specific material entropy production due to irreversible processes 
involving the climatic fluid, and T the temperature field. The following expression is usually 
adopted for s^at [14, 16, 27, 40]: 

Smat ='- + 'PsH-v(^+ ¥ lH " V [^ (2) 

where e^ the specific dissipation rate of kinetic energy, Ysh the turbulent sensible heat 
flux , Ylh the turbulent latent heat flux, where by turbulent we mean not related to large 
scale advection due to winds, which is in principle reversible. Romps [42] refers to the 
representation of the entropy production given by Eq. (2) as resulting from the bulk heating 
budget, because water is treated mainly as a passive substance, while processes such as 
irreversible mixing of water vapor are altogether ignored. More detailed description of the 
moist atmosphere have led to formulations of atmospheric thermodynamics able to account 
for these processes [8, 38, 39, 42]. Such more refined formulations of the entropy processes 
inside the climate system account for a consistent treatment of the entropy generated by the 
various processes accounting for hydrological cycle. Apparently, though, the overall effect 
of hydrological cycle-related entropy production is capture quite well using Eq. (2) [8, 29]. 

From Eq. (1), it is clear that there are two different ways to express the integrated 
material entropy production of the climate system, a direct one: 



^mat — I d ^ Smat (3) 

Jv 
where the specific material entropy production Smati which involves directly the properties of 
the climatic fluid, is integrated over the volume V , and an indirect one, involving radiative 
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heating rates only: 
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where obviously Sf^^^ = S^^t- Equation (1) provides an intimate link between the radiative 
fields and the material fiow properties inside the climate system. Moreover, Eq. (4) is very 
powerful because it permits to work out the material entropy production without passing 
through the details of the chemical and physical properties of the fiuid. The two equations 
(3) and (4) can be used to derive constants and bounds for linking the entropy production 
to other thermodynamical quantities of the climate system, as discussed in [29]. Pascale 
et al. [36] showed using the climate model FAMOUS adopted in this study that S^^f- and 
^mat agree up to an excellent degree of precision (within 1%). Since Pascale et al. [36] 
used the approximate expression (2) for the specific material entropy production, they also 
confirmed that, indeed, at all practical purposes using such simplified representation of the 
irreversibility associated to the hydrological cycle is appropriate. 



A. Entropy diagnostics in Climate Models 

In an actual climate model the implementation of entropy diagnostics faces some difficul- 
ties, both at theoretical level and in terms of practical implementation of the entropy-related 
diagnostics. 

A theoretical difficulty is that, as evidenced in [31], most state-of-the-art climate models 
features an inconsistent energetics, such that when all parameters are held fixed and the 
system reaches a steady state, the long-term average of the energy budget at the top of the 
atmosphere (TOA), which is the only boundary of the climate system, is unexpectedly biased 
with respect to the vanishing long-term average one should expect to observe. Interestingly, 
all biased models feature a positive energy budget at TOA, which implies that the time 
averaged outgoing long wave radiative fiux is smaller than the net incoming shortwave fiux. 
This fact implies that there must be a positive definite spurious sink of energy somewhere 
inside the system. More specific analyses make clear that such spurious sinks are related 
to the imperfect closure of the hydrological cycle [23] and to the inconsistent treatment of 
the dissipation of kinetic energy, which is not entirely (or at all) fed back into the system as 
thermal energy [31]. 
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These inconsistencies at small scales have relevant impacts on the large scale climatic 
properties: the positive bias in the TOA energy budget implies that the climate models are 
biased cold, taking into consideration that the Earth emits approximately as black body, 
or feature negative biases in the planetary albedo, or both. Moreover, since the biases are 
related to climate processes, they are climate-dependent, and so hard to control a posteriori 
via removal of anomalies. In terms of entropy production, an energy bias of the order 
oi 1 W ?n~^ causes a bias in the entropy production of about 4 x 10"'^ W m~^ K~^, 
which is comparable with the range of estimates of material entropy production even by 
various climate models [29, 36]. The FAMOUS model we use in this study features minor 
inconsistencies in terms of closure of the energy budget (the bias is smaller than 0.1 VT m~^, 
so that the problem exposed here does not affect significantly our results (see discussion 
later) . 

Moreover, in a climate model it is hard to deal directly with Eq. (2) because material 
turbulent fluxes are evaluated through parametrizations of unresolved processes. The cor- 
responding routines in the numerical code do not give as outputs heat fluxes . On the 
other hand the heating rates (i.e. the divergence of the heat fluxes) are easily diagnosed for 
these unresolved processes. Nonetheless, we can easily convince ourselves that using Gauss' 
theorem we can overcome this difficulty. A heating rate can be written in terms of the 
convergence of a flux (say e.g. F^h) as follows: 



over a volume V with boundary dV. If we integrate over the whole climate system, dV is the 
top-of-the-atmosphere and the bottom of the solid Earth. Neglecting the geothermal flux 
from the inner Earth and noting that at the top-of-the-atmosphere we have only radiative 
fields we conclude that the boundary term is zero. Hence, when the whole climate system 
is taken into account, the material entropy production can be worked out by considering all 
the material diabatic heating rates, as shown in [36]: 



3 I ^^\ (^'^Sh\ (^■Flh\_ f ,3_lQm.at 



St. - I d^x^- j - ^^^ j - ^^^) ^ I .'x^= j (6) 

In this paper we refer to the entropy budget of the FAMOUS atmosphere-ocean general 
circulation model [15] which has been studied in detail by [36, 37]. Lets first focus on the 
evaluation of S^^^. Different processes contribute to the entropy production terms described 
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in Eq.(2): the heating rates are calculated as output of many different parametrization 
routines describing the unresolved processes in the various subdomains of the climate system 
(atmosphere, ocean, soil, cryosphere): 



Entropy production due to dissipation of kinetic energy, Ske-, defined as: 



SKE=jd'^(^^y (7) 

In FAMOUS and, in general, in most climate models, the kinetic energy is dissipated 
mainly through four parametrized processes: the turbulent stresses occurring at the 
boundary layer, which extract kinetic energy from the free atmosphere, the gravity 
wave drag, which dissipates kinetic energy in the upper atmosphere, atmospheric con- 
vective processes, and small scale turbulence, which is represented by the horizontal 
momentum hyperdiffusion (which serves also the purpose of increasing the numerical 
stability of the model); only the atmosphere contributes to this part of the entropy 
production. 



Entropy production due to irreversible transfer of sensible and latent heat, Shf.at 



SsH + Slh, defined as: 



Shp.ni. — / c? X 



V-FsH\ / V-Flh 



- SsH + Slh (8) 

where we have separated the contributions coming from irreversible sensible and latent 
heat transfer. The boundary layer schemes contributes to the entropy production due 
to irreversible sensible and latent heat transfer in the four subdomains of the climate 
system, as it couples them through exchanges of sensible heat and water vapour; 



other processes contributing to Ssh and Slh are the atmospheric convection, whose 
parametrization describes the atmospheric vertical mixing in case of instability, the 
condensation and evaporation of water in the atmosphere, determined by the clouds 



and precipitation parametrization schemes, while processes contributing only to Ssh 
are the oceanic convection, the small scale turbulent mixing of temperature described 
by hyperdiffusion, and the mixing occurring inside the ocean associated to diffusion 
due to small scale eddies and due to the turbulent processes occurring in the mixed 
layer. One must note a basic difference between the contribution to entropy production 
due to sensible and latent heat irreversible transfer is that the former is eminently a 
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process where heat is absorbed and released in the same atmospheric column. In the 
case of latent heat, winds can transport away water vapor so that condensation of a 
water vapor parcel can take place at large horizontal distances from the region where 
evaporation had taken place. The closure of such mass transport through river and 
ocean transports is the essential element of the hydrological cycle. 

Table 1 provides a synthetic outline of which routines describing unresolved processes con- 
tribute to the various terms of the material entropy production in each climatic subdomain. 



Therefore, in practice, we compute S^l^^ as follows: 



•c 



k c -^^c ^ ^ 

where q1 is the local instantaneous heating rate occurring in the subdomain Vc due to the 
process k. 



The evaluation of S'^^at is much easier because the heating rates are readily available from 
the radiation scheme, which affects all the subdomains c of the climate system: 



. JVc 



Hsw \ I / Hlu 



T \T 



(10) 



where we have divided the contribution g^ui coming from the shortwave radiation, which 
is only absorbed (and scattered), inside the climate systems, so that q^w > 0, from the 
contribution g;^ coming from the longwave radiation, which instead is scattered, absorbed, 
and emitted, and is the sole responsible for the radiative cooling. 

III. COARSE-GRAINING OF THE ENTROPY PRODUCTION TERMS: DEFI- 
NITIONS AND SOME CONJECTURES 

As seen in the previous section, the entropy budget is estimated using space and time 
integrals of the ratio between the local heating term and the local temperature. In many 
cases, either because we need to compress data or because climatological database only 
contain certain time or spatially averaged data, we have to deal with coarse grained data 
for the heating rate (g(x, t))^, and for the temperature (r(x, t))^, where r refers to the time 
scale of the temporal averaging operation, and v refers to the set of stencil regions v{x) 

12 



centered over x over which (mass-weighted) spatial averaging is performed: 

{X{^,t))l = ^-— [d'y [ rfaX(x + y,t + a) (11) 

where fi{v{x),t) = j,sd^:K is the mass contained in the stencil v{x) at time t. Mass- 
weighting is the natural choice in climate models as hydrostatic approximation is almost 
invariably used and vertical coordinates are expressed to a very good approximation in terms 
of pressure levels. Since Eqs. (4) and (6) are not linear in the heating and temperature fields, 
we obviously have that for every r and v. 



cdir 


Jv \T J 
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// V (T)S )- 
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f ^^^f(<irad)l\ 

Jv \{T)l) 



Wmat/v^ \^^) 



Wmat/v- (13) 

Moreover, while as discussed before S^l^ = 5'™^^, there is no a priori reason to expect that 



{^'mat)l ^^d {S^^^)l have the same value. Finally, we have that up to first order: 
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cind /Cind\T _ A Q^'^d ^ / ^3^, [•^radj^ [ j^ /.. j.\ 

'-'mat Wmatlv ~ ^ '-'mat — / « X „ \^'^) 

L ^V Jy [(T)-] 

where A [Xj^ = X — {X)^. Equations (12)-(15) need to be considered when addressing 
the practical problem of estimating material entropy production from observational data 
or model outputs. They may address practical questions like: what is the error related to 
remapping the output of a climate model to a new resolution in space and time? How do di- 
urnal, seasonal and interannual variability and how different spatial structures (midlatitude 
cyclones, equator-pole contrasts, longitudinal asymmetries due to ocean-land contrasts, etc) 
affect the entropy budget? How should we proceed to compare the estimates of material 
entropy production from models with different resolutions? Note that this matter becomes 
extremely relevant when we consider the possibility of estimating the material entropy pro- 
duction for planetary systems where our ability of observing or modeling at high spatial 
and temporal resolution is pretty limited. These issues may be sequentially investigated by 
filtering g(x, t) and T(x, t) over the associated time- and space- scales. Moreover, we need to 
understand whether it is more accurate to obtain estimates of the material entropy produc- 
tion from coarse grained fields of the radiative heating rates or of the material heating rates, 
which can be used for the indirect or direct formula for the entropy production, respectively. 
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Given the sense of the operation of coarse graining, it is seems natural to interpret 



mat/ vi 



{^mat)l &s the contribution to the entropy production due to irreversible processes occurring 



qind 



A 



qdir 
'-'mat 



in Eqs. 



on scales large than what described by r and v. Consequently, A 
(14)- (15) can be interpreted as the contributions to the entropy production given by the 
material flows (Eq. (14)) and radiative fluxes (Eq. (15)) with variability confined below the 
spatial scale given by v and by the time scale given by r. 

When we perform the coarse graining given in Eq. (11) to the thermodynamic variables 
and estimate the entropy production, we discount for the mixing processes occurring below 



the chosen spatial and time scales. Therefore, one expects that (S'^aJjJ; {Smat)li ^ ^^^ 

> for all choices of r and v. Moreover, it seems natural to conjecture 



A 



qind 



,A 



qdu- 



that if, given a model output, we choose a coarser graining, we should obtain a lower 
estimate of the entropy production, because we neglect a larger set of irreversible processes. 



In other terms, we should have that A 



qdir 
'-'mat 



> A 



>b,„ 



mat 



[or (Omai)t;l < \Srnat)v\) and 



A 



qind 
'-'mat 



> A 



qind 



mat 



J V2 



(or {St^i)l\ < {St^it)Z) if ^2 < n and V2 ^ v^. Let's see how to 



interpret the inequalities A 
• The inequality A 



qind 
'-'mat 



A 



qdir 
'^mat 



> using the r.h.s. of Eqs. (14)-(15): 



qind 
'-'mat 



> can be interpreted as the fact that at all time and space 
scales, there is on the global average a positive correlation between the anomalies of 
radiative heating and the anomalies of temperature. This expresses the basic fact 
that the climate system is driven by radiative forcings, in the first place, and this 
terms refers to the response of the system to the external forcing. Note that the 
correlation between the absorbed shortwave and longwave radiation and the temporal 
field must overcome the strong negative correlation due to the fact that any body 
emits proportionally to the fourth power of its temperature. 



The other inequality A 



qdir 
'-'mat 



> 0, instead, implies that at all time and space scales 
on the average there is a negative correlation between the anomalies of heating due 
to convergence of material heat fluxes and anomalies of temperatures. This relation, 
instead, expresses the fact that temperature anomalies are damped by the geophysical 
flows, and this terms refers to the dissipation occurring inside the system. 

In other terms, these inequalities clarify the well-known fact that the climate is forced 
by anomalies in the radiative forcing and the atmospheric and oceanic circulation 
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reduce the resulting temperature gradients. 
In this perspective, it reasonable to expect that 

A 



Qind 



mat 



> A 



Cdid 
'-'mat 



Vr,w, 



because otherwise the system would be able to dampen perfectly all external forcings. There- 



fore, we can conjecture, considering Eq. (14)-(15) and the fact that 5"^^^ = S^^^ti that 



\*-'^^/i//^, ^ \ *-'^>i/-,+ /io vT.V 



\'~^matlv — \'~^mat/v 



A. Some examples 

Before proceeding further, let us consider some limiting cases for Eqs. (12)-(13). Let 
us consider as averaging volume v a vertical column ranging from the bottom of the fluid 
component of the climate system to the top of the atmosphere, and let us consider a long 
averaging time r > 1 y, so that all temporal dependencies are removed. 

We start with the expression relevant for the indirect formula for estimating the material 
entropy production: 



{<irad)l\ f, , Ftoa(xi,X2) 



(«.); = / d^x( ^:^ ) = - / dxidx.l^ll^^il^i^. (16) 

where Ft oAi^i 1^2)1 resulting from the application of Gauss' theorem, is the climatolog- 
ical average of the net radiative wave flux at the top of the atmosphere (positive when 
there is net incoming radiation towards the planet) at the horizontal coordinates Xi,X2 
, while Tc/j(xi,X2) is the long term mean of the vertical average of the fluid tempera- 
ture. Such quantity can be closely approximated by the emission temperature T£;(xi, X2) = 
(F^o^(xi, X2)/cr)-'^/'^, where a is the Boltzmann's constant and F^q^(xi,X2) is the outgoing 
long wave radiation at TOA. We remind that Ftoa(xi, X2) = F^q^(xi, X2) — F^q^(xi, X2), 
where 

Note that since Ftoa and Te are positively correlated (regions having a net positive 
incoming radiation are warmer), we have that Ftoa and l/T^ are negatively correlated. 



Since j^ dxidx2-FroA = 0, we derive that {S'!i^at)l i^i Eq. (16) is positive, as expected. 
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Equation (16) can also be given a different interpretation. We have that Fto^(xi,X2) 
is equal to the divergence of enthalpy transport due to the large scale climatological atmo- 
spheric and oceanic flow, so that 



i-TOA 



/lUA 
dx3[YL/^(x) + Tsh(x) + Tpot(x)] 
- sur f 

= V2 ■ [Tlh{xi,X2) + Tsh{xi, X2) + Tpot{xi, X2)]. (17) 

In the previous equation, we have indicated with Ylh{^) = L^qi^x) {vi{'x) , V2{x.)} and 
T5/^(x) = CpT(x.) {t;i(x),f2(x)} the large scale, advective horizontal fluxes of latent and 
sensible heat, respectively, where L^ is the latent heat of evaporation of water (taken as 
a constant for simplicity), q is the specific humidity, and Cp is the heat capacity of air at 
constant pressure, and Vi and V2 indicate the two components of the horizontal velocity 
field. Instead Tpot(x) = gx2,{vi{yi),V2{'x)} represents the large scale horizontal advection 
of potential energy [31, 40], is much smaller than the other terms and can be neglected 
[40] . Finally, the ~ sign refers to the vertically integrated fluxes. Inserting the right hand 
side of Eq. (17) in Eq. (16), we conclude that Eq. (16) can be interpreted as the entropy 
produced by the large scale horizontal transport of the geophysical flows and describes the 
climate system as a purely 2D system featuring irreversible heat transport from warm to 
cold regions. 

Let's now repeat the same coarse graining operation for the direct formula. As discussed 
above, in each location, when integrating vertically, the surface sensible heat fluxes cancel 
out with the heating rates associated to sensible heat fluxes in the atmosphere, so that for 
this choice of coarse graining their contribution to the entropy production vanishes. We 
obtain: 



{S^'at)l = I d^X^i^:^) = 1^ dXidx, 



L^,[P{xi,X2) - E{xi,X2)] e2(xi,X2) 



dxidx2 



rdi(Xl,X2) Tcii(Xi,X2) 

V2 ■ tL/^(xi,X2) , e2(xi,X2) 



^di(Xl,X2) Tc;j(Xi,X2j 



(18) 



where Lw[P(xi,X2) — E(xi,X2)] = — V2 ■ Ti//(xi,X2) as imposed by conservation of water 
mass [40], we indicate with e^(xi, X2) the vertically integrated kinetic energy dissipation rate, 
and we choose, as a first approximation Tc/j(xi,X2) as characteristic temperature defined as 
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before. Taking the difference between Eq. (18) and Eq. (16), we obtain: 

e^(xi,X2) , +V2- t5j/(Xi,X2) 



{SiZatY. - {St:ii)l ~ j dXidX2 



+ 



T'cZi(xi,X2) Tc;i(xi,X2) 



> 0. (19) 



where the inequahty, which agrees with our previous general conjecture, derives from what 
presented in [29], where it is proved that the entropy produced by dissipation of kinetic 
energy is larger than what produced by large scale material transport of latent and sensible 
heat (here only sensible heat comes into play). Our previous conjecture suggests that this 
quantity is positive. This is actually the case, as Moreover, in [29] it was derived that for 



this choice of coarse graining {S!^at)l — {^matY ^^ ^^Y climate, so that the coarse grained 
direct and indirect formulas are not equivalent. Specifically, using parameters suitable for 
the present terrestrial climate, we get that Eq. (16) accounts for about 15% of the total 
material entropy production on our planet, while Eq. (18) accounts for about twice than 
that [29]. 

We can bring the previous example to a more extreme case. If by v we consider the whole 
climate domain, it is easy to derive that: 



V^— , ^ {qrad)l , ^ /sdXidX2FTOA(xi,X2) 
mat)l=V = ^(^)^fyr = -^(^) ^fp = °' (2°) 

because through such an extreme coarse graining we have reduced the climate to a zero- 
dimensional system with a unique temperature where absorbed and emitted radiation are 
equal. Such a system is at equilibrium, cannot do any work, and cannot sustain any irre- 
versible process. 



Interestingly, if we compute {S'^at)l using the same global coarse graining, we do not 
obtain a vanishing result. While the contribution to entropy production due to heat fluxes 
is eliminated, because we integrate over a domain the divergence of fluxes which vanish at 
the border, the contribution coming from the dissipation of kinetic energy is 



{qmat)l _ , ./sdXidX2?(xi,X2) 



{S^:t)l=v = f^i^')^==r = M^) ^" ' ^' '' " > 0. (21) 

Equation (21) gives, to a good degree of approximation, the minimum value of the entropy 



production compatible with the presence of a total dissipation jy (i^xe2(x, t) [27, 29]. There- 
fore, the physical interpretation of Eq. (21)is rather different from what given in Eq. (20), 
even if the coarse graining is the same. 
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IV. RESULTS 

We first discuss briefly how the coarse graining operation is performed in practice. Let 
us consider a steady-state chmate simulation lasting for a time period L (in our case L = 50 
years), which we divide it in A^ sub-intervals r = L/N , where r = M x dt, where dt is the 
model's time step (1 h in our case). The horizontal resolution is specified by regular grids 
with angle resolution of 5° x 7.5° lat-lon), while in the vertical we have 11 levels for the 
atmosphere, 20 oceanic levels, and 3 land surface levels [15]. Therefore, we subdivide the 
domain V of integration into Q subdomains Wg, g = 1, . . . , Q, each containing (in the bulk of 
the model's domain) R grid points. Given an intensive thermodynamic field X{xk,tj), for 
n = 1, . . . , A^ and q = 1, . . . ,Q, we define its coarse grained version as: 

-. Rq Mn 

(X(g,n)); = -— — Y. E dtMx„a,)X(x,,a,), (22) 

'^^ •?' ' j=ij(q_i)+i k=M{n-l)+l 

where /u(xj, cr^) = i/(xj)p(xj, 0"^) is the mass contained in the grid box centered around Xj of 
volume i^(xj) at time at and /u(wg, n), correspondingly, is the time averaged (for time ranging 
from tM(n-i)+i and tun ) niass contained in the domain Vq of volume vivq). Therefore, our 
estimate of the coarse grained value of the material entropy production is: 

^S^^=^EE-(-^)%f^ (23) 

for the so-called direct formula, and: 

for the so-called indirect formula. These formulas are the discrete versions, implemented as 
they are in the post-processing of the FAMOUS climate model, of the continuous equations 
given in the r.h.s. of Eq. (12) and Eq. (13), respectively. Obviously, the discrete versions 



of the exact formulas for the material entropy production 5**^^ and S^^^ are obtained by 
setting in Eq. (22) i? = M = 1, i.e., taking the model outputs at the highest possible 
resolution. 

If we choose R = 1 and we consider different values of M, we test how applying temporal 
coarse graining only impacts the estimates of the material entropy production. Instead, if 
we change the shape of the stencil v and/or the number of points R is contains while keeping 
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M = 1, we study how different spatial coarse graining alter the estimates of the material 
entropy production. We have to remember that we cannot capture the contributions to the 
material entropy production due to irreversible processes taking place over timescale shorter 
than the model timestep and over space scales smaller than the model resolution. It is not 
clear, given a specific model's settings, how relevant these could be, and, indeed, the only 
way to find this out is to alter the model's resolution. This procedure may have relevance in 
terms of model tuning, as one could decide to change a model's parameter when altering its 
resolution in such a way to keep the entropy production constant. We also want to remark 
that we exclude from the computation of the entropy production the processes occurring in 
the interior of the ocean and below the first soil level. The corresponding terms have been 
shown to be entirely negligible [37]. 

A. Temporal coarse graining 

We start our investigation by performing coarse graining exclusively on time. We then 
analyze a long, steady state model's run lasting 50 years with a model's timestep of 1 hour, 
and consider 1 year as long-term averaging time. We use the following values for r: 1 hour 
(model timestep, iV = 1), 6 hours (A^ = 6) , 12 hours (A^ = 12), 1 day (A^ = 24), 2 days 
(A^ = 48), 5 days (A^ = 120), 10 days (A^ = 240), 15 (A^ = 360) days, 1 month (A^ = 720), 
3 months (A^ = 2160), 6 months (A^ = 4320), 1 year (A^ = 8640). We the collect the 50 
1-year averaged value of the coarse grained material entropy production and compute the 
mean and standard deviation for the 50 data we have. These data are collected to analyze 
the properties of the climate system for scales ranging up to 1 year. Moreover, we consider 
longer averaging periods - 5 years, 10 years, and 50 years, and take in these cases r equal to 
the averaging time, so that A^ = 43200, A^ = 86400, and A^ = 432000 in the r = 5, 10, and 
50 years case, respectively, thus spanning in total more than 5 orders of magnitude for A^. 
We then compute for the coarse-grained estimates of the material entropy production the 
ten 5— year averages and the five 10— year averages, and compute the mean and standard 
deviation, plus the unique value referred to the 50— year average. 

In Fig. 1(a) we report the estimates of the material entropy production obtained through 
the direct formula {Sf^l^^Y and the indirect formula {Sf^^^y , respectively, where we have 
dropped the lower index v because we do not perform any spatial coarse graining. The 
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vertical bars indicate the uncertainty due to the long term variability. 



The computed values (worked out at each timestep) of S'^fj f« 53.1 mW m ^ K ^ and 



^mat ~ 53.5 mW m~^ K~^ have a difference of about 0.4 mW m~^ K~^, so that Eq. (1) is 
verified with great accuracy. The discrepancy term between the two estimates is due to the 
extremely small spurious radiative imbalance at TOA of about 0.1 VT m~'^ (see [29]) and 
to numerical inaccuracies. Moreover, as discussed in [29, 37], these estimates are in good 
agreement with what found in climate models of higher degree of complexity. 

As expected, we find that the estimates of the coarse grained entropy production decrease 
with increasing r from these reference values obtained with no temporal coarse graining. As 
conjectured, the bias resulting from the use of the indirect formula is larger for all values of 
r. In Fig. 1(a) we see that if we consider values of r up to 6 hours, the impact of coarse 
graining is negligible. This implies that such small time scales the irreversible processes 
are negligible; this matches well with the fact that convection, which is the dominating 
fast process in the climate system, is parametrized with an instantaneous adjustment. This 
immediately points to an unwelcome spurious effects of climate parametrizations. 

The effect of coarse graining becomes more relevant when r > 1 day. Figures 1 (b) and 
1(c) present the values of A 



Qdir 
'-'mat 



and A 



Qind 



as a function of r. As we see, for r 



■^mat 

day, the difference between the true and the coarse grained value of the entropy production 
is about K. 0.6 mW m~^ K~^ if we use the direct formula, and ~ 2 mW m~^ K~^ is we use 
the indirect formula. Such biases are due to neglecting the mixing occurring on the time 
scale of the day, and related to the daily cycle of incoming radiation. When considering the 



direct formula, it is interesting to note that A 5*^'^ is basically zero for all values of r (not 
shown), meaning that there is no time correlation between the dissipation of kinetic energy 
and the temperature field. The dissipation of kinetic energy is usually small, and does not 
influence strongly the local temperature; nonetheless, it plays an important role because, as 
discussed above, it is a positive definite contribution to the entropy production which resists 
all coarse graining procedures. What instead we are neglecting is the effect of the diurnal 
cycle on the hydrological cycle (e.g. via convection). We can substantiate this statement by 

-| T r -| T 

observing that A S^q^ ~ A 5'^g'Jjj , as can be seen by comparing Figs. 1(c) and 2(a). 



The second timescale worth discussing is the one corresponding to 1 year (~ 3 x 10^ s). 
The use of annual means instead of time-step data introduces a further bias of about 2 mW 
m~^ K~^ when using the indirect formula, which corresponds to neglecting the correlation 
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between the seasonal cycle of the radiation budget and that of the radiation temperature 
field. Similarly, considering the direct formula, an additional bias of about 1.5 mW m~^ K~^ 



given by A 5*^^ estimates the effect of neglecting the correlation of the seasonal cycle of 
the atmospheric and oceanic transport and dissipation, and that of the temperature field. 

for r > 1 year 



One must note that a considerable contribution to the value of A S'^^j 
is given by the atmospheric temperature hyper diffusion, which in FAMOUS is implemented 
as a eight-order laplacian operator and applied after the advection to the model prognostic 
variables. Hyper diffusion is generally introduced in dynamic cores for numerical reasons in 
order to smooth variables and avoid local divergences. However it may thought as a way to 
represent turbulent dissipation and mixing at subgrid scale. Also in this case, we discover 
that a traditional numerical trick used in the climate modeling community for avoiding 
computational troubles has a considerable large scale impact on the budget of physical 
quantities of the system. 

We also observe that there is no clear signature emerging in the functions A 
A 



qdir 



Cind 

^ Til ni 



and 

for values of r to weeks (synoptic waves) or monthly (low frequency variability) 



time scales while a relatively smooth transitions is found going from daily to yearly averages. 
This supports the idea that it is possible to look at weather disturbances as parts of a macro- 
turbulent cascade. 

Therefore, our analysis quantifies the impact of the two main time modulations of the 
solar forcing on the climate system on the entropy production. We discover that estimating 
the entropy production via either the direct or the indirect formula using long term averages 
(but full spatial resolution) we underestimate the exact value of the entropy production by 
less than 10%. This suggests that long term averages of the climatic fields one can obtain 
from the climate repositories are enough to get a good idea of the properties of the climate 
system. It is crucial, though, to keep the information on the spatial structure of the heating 
and temperature fields. 

We conclude this section with a note on the oceanic processes, which we do not treat 
in this paper as they contribute negligibly to the overall entropy production in the climate 
system, in Fig. 2(a) we show the dependence of the entropy production due to the oceanic 
mixing on the temporal coarse graining (dashed) line. We discover that its exact value, 
computed at time step, is about 1 mW m~^ K~^, as in [37], and its coarse grained value 
does not noticeably decrease up to r ~ 1 year, above which the coarse grained estimate is 
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roughly halved. The dash-dotted line in Fig. 2(a) gives the contribution due to the vertical 
mixing in the interior of the ocean, which is a very slow process and is, in fact, weakly 
affected by the temporal coarse graining. The other contribution to the entropy production 
in the ocean comes from the mixing occurring in the mixed layer. The mixing layer scheme 
[20] parametrizes the convection due to heating at depth and cooling at the surface as well 
as the mechanical stirring due to wind and is introduced in ocean models in order to account 
for the seasonal thermocline variations. The coarse grained value of this term goes virtually 
to zero for r > 1 because, when considering such an averaging, we discount for the impact 
of the seasonal cycle in the upper portion of the ocean. 

B. Spatial and Temporal coarse-graining 

In this section we analyze the combined effect of coarse graining the heating rates and 
temperature fields in space and time by using extensively Eqs. 23 and 24. We remind that at 
full resolution the spatial discretization of the atmospheric variables is described by a regular 
grid with angle resolution of 5° x 7.5° lat-lon and by 11 vertical levels is a coordinates. Of 
course, there are many ways to perform coarse graining, boiling down to the selection of the 
stencil v introduced before. Summarizing, we proceed as follows with the coarse graining: 

1. longitudinal averaging: the stencils v are given by arcs of varying length in the zonal 
direction; 

2. areal averaging: the stencils v are given by portions of varying size of the spherical 
surface; 

3. mass averaging: the stencils v are given by same-mass portions of the atmospheric 
spherical shell obtained by thickening in the vertical direction the stencils described 
in 2.; 

in all cases we perform also temporal coarse graining by selecting the same averaging times 
r described in the previous subsection. We underline that, given the resolution of the model 
output, choosing stencils v with the same mass throughout the volume of the climate system 
would require heavy remapping on the thermodynamic fields, which would lead to spurious 
contributions to the entropy production. Instead, we choose coarse grained grid boxes that 
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respect as much as possible the original model's resolution. We also remark that given the 
heavy computational burden of the operation, we restrict our analysis to only one of the 
fifty year of available data. 

1. Longitudinal Averaging 

We first investigate the effect of coarse graining on the estimate of the material entropy 
production by averaging longitudinally the thermodynamic fields, up to the point of con- 
sidering zonally averaged only fields, and by degrading their temporal resolution by using 



the averaging times r described above. The estimates of (S^^Jj; and {S^at)l ^^^ given in 
Fig. 3(a) and Fig. 3(c), respectively. In both figures, the value of r is reported in the abscis- 
sae, while in the ordinates the value of size of the spatial stencil, ranging from 7.5° (no coarse 
graining) to 360° (zonal averaging) is shown. In both figures, the lower left corner corre- 
sponds to the best estimate of the entropy production; the upper left corner corresponds to 
the material entropy production due to the longitudinally averaged, high-temporal resolution 
fields, the lower right corner corresponds to long-time averaged, high-resolution spatial case, 
and, eventually, the upper right corner corresponds to the highest degree of coarse graining: 
it represent the entropy production due to the long-term averaged, longitudinally averaged 



fields, and features the lowest value of {S'^^t)l and {S!^at)l- We remark that the values 
reported at the border of the domain given by the lowest value of the ordinates coincide, ob- 
viously, with what shown in Fig. 1(a). As a general fact, we observe that the coarse grained 
estimates of the entropy production decrease (or remain virtually unchanged) as we perform 



coarser and coarser graining procedure, in time or in space, and that {S^lt)v — {^mat)l- 



The strongest dependence of the {S^at)l i^ on r: temporal coarse graining appears to 
be the dominating influence, while the effect of spatial coarse graining is apparent only 
for r < 1 day and for considerable longitudinal averaging, such that features below 60° 
are smeared out. In other terms, longitudinal averaging starts to matter only when we 
lose information on the alternating pattern continents/oceans. The total effect of removing 
totally the spatial structure is similar to that of performing a time-averaging of one day. 



When considering {S^a^)l the picture is partially different: first, the influence of the spatial 
averaging is relatively strong at all scales for r < 1 day. The coupling between the spatial 
and temporal scales indicates that using low pass filter and space and time we remove the 
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fast traveling synoptic waves of the mid-latitudes. As opposed to the case of the coarse 
grained indirect estimate of the entropy production, spatial averaging plays a role also for 1 
day < r < 3 months. This is probably the signature of the relevance of low-frequency, large 
scale features of the tropical circulation, which are sustained by longitudinal gradients (and 
tend to reduce them), which are smeared out when extreme coarse graining is applied. 

Concluding, we remark that neglecting information on the longitudinal fluctuations and 
temporal fluctuations of the thermodynamic fields does not bias considerably (in the worst 
case, by about 10%) the estimates of the entropy production one would obtain by retaining 
the full information. This agrees with the fact that, in first approximation, in our planet 
longitudinal gradients and longitudinal heat fluxes are relatively small [40] . Moreover, using 
either the direct or indirect formula we obtain rather similar results, with a bias of maximum 
5%. As expected, the direct formula gives more accurate estimates for all considered coarse 
graining procedures. 

2. Areal averaging 

As a second step for understanding the role of spatial and temporal coarse graining of 
the estimate of the material entropy production, we combine time averaging of the ther- 
modynamic fields with areal averaging along horizontal surfaces. This operation allows us 
to explore how the two-dimensional spatial covariance of heating and temperature fields 
contributes to entropy production at all time scales. In order to keep coherence with the 
previous coarse graining procedure, we proceed as follows. We divide the spherical surface in 
coarse grained grid boxes defined by intervals (in degrees in latitude and longitude) (AAA0) 
such that A0/AA = 1.5, which is consistent with the model's resolution of 5°lat x 7.5°lon. 
We then increase A0 from 7.5° up to 90°, thus decreasing progressively the number of coarse 
grained grids from 1728 to 12. In order to complete the coarse graining, we select as last 
two resolutions (AA; A0) = (90°, 180°) (four quadrants) and (AA; A0) = (180°, 360°) (full 
sphere). 



The estimates of {S^'^f)l and {S^at)li ^^^ given in Fig. 4(a) and Fig. 4(c), respectively. We 
discover that, as opposed to the previous case, the impact of selecting coarser and coarser 



graining in space reduces considerably the value of {S^ai)l for all values of r, because such 
an averaging progressively removes the strong meridional dependence of the thermodynamic 
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fields, up to the extreme case of v being tlie wliole Eartli's surface. In this case, the estimate 



of the entropy production {Sflll^)l ~ 47.2 'mWm~'^K~^. Note that when considering very 
strong spatial averaging the effect of changing r is negligible, because the spatial averaging 
alone reduces the temporal correlations by mixing areas of the planet experiencing, e.g., 
different seasons. The r-dependence of {S^lt)l is relevant only for r < 1 day and spatial 
scales smaller than ~ 3 — 4 x 10^ m, which is, like in the previous case, hints at the 
fact that when averaging over large spatial and temporal scales, we remove the variability 
corresponding to synoptic waves. 



The function {S^mat)l has a qualitatively similar but quantitatively stronger dependence 
on A0 and r with respect to {S^tt)l'- ^^e bias in the estimate production is larger for all 
the considered coarse graining. Additionally, the indirect formula is more strongly affected 
by of averaging over long time scales r, similar to what seen in Fig. 3(c), because the 
coupling between the seasonal cycle of the radiative budget and the temperature fields is 
very strong. Note that when we consider global or quasi-global spatial coarse graining, such 
effect disappears, because averaging such large scales already removes a large part of the 
season cycle signal. This is different from what reported in Fig. 3(c), because zonal averaging, 
obviously, cannot remove the asymmetry between northern and southern hemisphere. 

3. Mass Averaging 

Finally, we perform spatial averaging along the horizontal and vertical direction and for 
different values of r for the direct and indirect formula of entropy production. In this way, 
we are able to ascertain the relevance of the processes involving irreversible fluxes across 
temperature gradients along the vertical direction. Results are reported in Fig. 5 and Fig. 
6, respectively. Since we are now dealing with three variables describing the coarse graining 
- the amplitude in latitude Ac/), the number of levels n, and r, we present two cross sections 
obtained for r = 1 day (virtually indistinguishable from r = 1 hour) and r = 1 year, 
reported as panels (a) and (c), respectively, in both Figs. 5 and 6. In panels (b) and (d) of 
these two figures, we report, instead, A 



qdir 



and A 



Qind 



respectively. The inequality 



A 



Cdii 



mat 



< A 



J V 



Cind 



mat 



is clearly obeyed. 



J V 



As we know from the previous discussions, the function {S^'!^f)l is weakly affected by 
coarse graining along the horizontal directions and along the time axis. The modest im- 
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portance of time averaging is confirmed in this more complete analysis, as Figs. 5(a), and 
5(c) are hard to distinguish. We find that averaging the thermodynamic fields along the 
vertical reduces the temperature gradients across which the irreversible heat fluxes transport 
heat, thus reducing drastically estimate of the entropy production, in a quasi-linear fashion 
with respect to the number of levels over which the averages are performed. The effect of 
horizontal averaging becomes relevant only when we increase the effective A0 to 90°, so that 
the horizontal gradients are almost entirely smeared out. 

When considering the indirect formula, we can draw roughly the same conclusions as 



above, with the difference that {Sl^^i)l i^ much more strongly affected by averaging along 
the horizontal surface and along the time axis. Therefore, Figs. 6(a) and 6(c) feature clear 
differences in terms of mean values, and, in each of them, the impact of performing very 
coarse graining along the horizontal direction is more pronounced with respect to what 
reported in Fig. 5. Moreover, the effect of horizontal coarse graining become noticeable 
already when grid boxes with side of (A A, A0) ~ (20°, 30°) are considered. 

4-. Comments 

As previously discussed, performing coarse graining to the estimate of the entropy pro- 
duction can be interpreted heuristically as neglecting irreversible processes occurring below 
the scales over which averaging is performed. One has to keep in mind that, nonetheless, 
when considering the direct formula for the estimation of the entropy production, the contri- 
bution to due to the dissipation of kinetic energy is virtually unaffected by coarse graining. 
Let' se see the conclusions we can draw from this point of view. We consider Tables II 
and III, where the values of A 



qdir 
'-'mat 



and A 

J V 



qind 
'-'mat 



obtained for extreme coarse graining 



procedures are reported. 



Along these lines, we can then interpret A 5*^^^ = 5.3 mWm ^K ^ in Table II {Areal 



TM 



VM 



averaging row) as the portion of entropy production due to climatological sensible and latent 
heat transport along horizontal surfaces. Looking at the row above [Longitudinal averaging 
= 2.2 mWm~'^K~^ as the the entropy produced by the climatological 



row), A 



qdir 
'-^mat 



horizontal transport along the zonal direction, i.e. across the land-sea temperature contrast 
. The difference between these two values can then be seen as the contribution due to 
climatological meridional heat transport (= 3.1 mWm~'^K~^). Somewhat surprisingly, the 
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contribution of the meridional heat transport is only 1.5 times that of the zonal transport. 

We can check this interpretation using a different approach. When we consider large val- 
ues of r ~ 1 year and vertical averaging along the full atmospheric column, while performing 



no coarse graining in the horizontal direction, {S'^mit)v corresponds to the entropy produced 
by the irreversible heat transport due to the climatological large scale flows, as discussed in 
Eq. (16). Looking in the bottom right corner of Fig. 6(c) (which corresponds to the last 
entry in Table III), we obtain a value of ~ 6.5 mWim~'^K~^, which is reasonably close to 
the value of 5.3 mWm~'^K~^ mentioned before. 

Further information can be derived from these figures. Equation (18) suggests that 
the bottom right corner of Fig. 5(c) should provide a value corresponding to the sum of 
entropy produced by large scale transport of latent heat plus the entropy produced by the 
dissipation of kinetic energy. We find a value of ~ 18 mWm^'^K"^. Considering that the 
entropy production due to large scale transport of sensible heat is much smaller than the 
corresponding contribution due to latent heat transport, we can derive that the dissipation 
of kinetic energy contributes for about ~ 13 mWm~'^K~^ to the total material entropy. This 
value agrees very well with what obtain by direct estimate of Ske ~ 13.5 'mWm~'^K~^. 

As additional test of consistency, we have that considering the coarsest graining possible 
in our model (f=full climate system, r = 1 year), we obtain, consistently with Eq. (20), 



that {S'^matll ~ (see Fig- 6(c) and penultimate entry in Table III). 

These results imply that we can obtain an extremely good estimate of the true value of 
the material entropy production even using climatological, horizontally global averages of 
the thermodynamics fields: in such a worst case scenario, we get a bias of about 10%. Using 
few selected coarse grained estimates for the entropy production, one can derive that it is 
possible to split the total value of 52.5 mWm~'^K~^ as follows: ~ 6 mWm~'^K~^ can be 
attributed to large scale transports of latent and sensible heat (of which ~ 2 mWm^'^K"^ 
due to zonal transports); ~ 13 mWm~'^K~^ can be attributed to the entropy produced by 
dissipation of kinetic energy; ~ 33.5 ■mWm~^K~^ can be attributed to vertical transports 
of sensible and latent heat (basically, convection). These estimates agree quite accurately 
with what obtained in [29] using scaling analysis. Also, one obtains, in agreement with the 
inequality proposed in [29], that the entropy produced by dissipation of kinetic energy is 
larger than that due to large scale heat transport. 
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V. CONCLUSIONS 

The investigation of the chmate system using tools borrowed from non-equihbrium ther- 
modynamics [2, 19] is a very active interdisciphnary research, which allows for connecting 
concepts of great relevance for climate dynamics, such as large scale heat transports and 
the Lorenz energy cycle [41], to basic thermodynamical concepts used in the investigation 
of general non-equilibrium systems [16, 27]. Such a theoretical framework seems relevant 
especially in the context of the growing field focusing on the study of the atmospheres of 
exoplanets, for which detailed measurements are hardly available. In fact, thermodynam- 
ical methods allow for defining inequalities and deducing apparently unexpected relations 
between different physical quantities [29]. 

In previous studies, we have emphasized how the material entropy production is an ex- 
tremely useful thermodynamical indicator of the state of a planet system, and provides, e.g., 
the clearest indication of whether a planet is a snowball state [1, 30]. We have also suggested 
, performing scale analysis over a formula provided by Goody [8], that it is possible to es- 
timate the total material entropy production and its portion due to large scale horizontal 
heat transports using low resolution data, and proved the effectiveness of this approach on 
actual climate models [29]. 

In this paper we have focused on understanding how the different spatial and temporal 
scales of motions contribute to the material entropy production in the climate system. We 
have considered the output coming from a 50 y run under steady state conditions performed 
with the FAMOUS climate model [15], and have used the entropy diagnostics developed and 
tested in [36]. We have considered both the direct and the indirect formulas for material 
entropy production [8]: the former estimates the material entropy production using the 
heating rates associated to the dissipation of kinetic energy and the convergence of material 
heat fluxes, the latter uses, instead, the heating rates associated to radiative fluxes. 

Our strategy has been the following: we have considered the estimates of the entropy 
production coming from the coarse grained outputs of the heating and temperature fields. 
The coarse-graining has been performed both in time and in space. The temporal coarse 
graining ranges from hourly (timestep of the model) to yearly time scale, while the spatial 
coarse grained has been performed in three different modalities: 1) performing longitudinal 
averages; 2) performing averages along horizontal surfaces; 3) performing mass-weighted 
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averages along the horizontal and vertical directions. We have discovered that the coarser the 
graining of the data, the lower is the resulting estimate of the material entropy production, 
both in the case of the direct and of the indirect formula for estimating the material entropy 
production. This implies that at all scales there is a negative correlation between heating 
rates related to flow (kinetic energy dissipation, sensible and latent heat fluxes) and the 
temperature fleld, and a positive correlation between the heating rate due to radiation and 
the temperature fleld. In other terms, at all scales, the climate systems results to be forced 
by radiation, while the resulting forced fluctuations are dissipated by the material fluxes. 
In agreement with this interpretation, at all scales the correlation between the radiative 
heating and the temperature fleld is stronger than the correlation between the temperature 
fleld and the material heating rates. If the two correlation were equal, the climate system 
would be able to adjust, instantly and locally, to (spatial and temporal) variations in the 
radiative heating. 

Considering various special cases of coarse graining, and using the basic thermodynamic 
equations, we have been able to estimate in a consistent way the contributions to material 
entropy productions coming from large scale horizontal heat transport (~ 6 mWm'^^K"^), 
dissipation of kinetic energy (~ 13 mWm?~K~^), and vertical processes of sensible and 
latent heat exchanges (i.e. convection, ~ 33.5 mWim?~K~^). This suggest that, as flrst 
approximation, the climate system can be seen in terms of dissipative processes as a collection 
of weakly coupled vertical columns featuring turbulent exchanges and dissipation. 

Note that one could use the quantitative information on the various contributions to the 
material entropy production to derive some basic properties of the climate system without 
resorting to the full three dimensional, time dependent flelds. In particular, one can derive a 
good estimate of the intensity of the Lorenz energy cycle by multiplying the estimated value 
of the contribution of the dissipation to the kinetic energy to the material entropy production 
times a characteristic temperature of the system, obtaining an estimate of ~ 3 Wm~'^, which 
is in good agreement with what obtained after processing of the high- resolution data [36] . 

The fact that the estimate of the material entropy production of the climate system 
decreases when a coarser gaining is considered is in qualitative agreement with what derived 
in the Appendix A for the simple case of continuum systems featuring generalized flux- 
gradient relations. This obviously does not imply that the climate system behaves as a 
diffusive system, yet share with a diffusive system this interesting property. The fact that at 
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all spatial and temporal scales the system has a positive definite value of material entropy 
production, when global averages are considered. This is not in contradiction with the well- 
known phenomena of so-called negative diffusion, first noted by Starr [44], who observed 
that in certain portions of the atmosphere - namely, near the storm track - the fiuxes of 
momentum transport momentum from low to high momentum regions. While this process 
- a crucial element of the general circulation of the atmosphere, observed in our model 
runs as well - seems to oppose the second law of thermodynamics, it is instead a local but 
macroscopic phenomenon, where the creation of organized structures (thanks to long-range 
correlations due to wave propagation), is, as we understand in this paper, over-compensated 
by large entropy production at the same scales elsewhere in the atmosphere. 

Apart from providing insights on the properties of forced fiuctuations and irreversible 
dissipative processes in the climate system at various spatial and temporal scales, this paper 
serves the purpose of providing information on what is the impact of being able to access 
data at lower resolution with respect to the model which has generated them. We have 
learnt that this lack of information always biases negatively our estimate of the entropy 
production, and that the bias is serious only if we miss information describing the vertical 
structure of thermodynamic fields. 

Since performing time averages up to the yearly time scale does not bias substantially 
the estimates of the material entropy production, we have that it is possible to intercom- 
pare robustly the state-of-the-art climate models and assess on each of them the impact of 
climate change on the entropy production by resorting to the output data provided in the 
freely accessible PCMD1/CMIP3 (http://www-pcmdi.llnl.gov/ipcc/about_ipcc.php) 
and PCMDI/CM1P5 (http://cniip-pcnidi.llnl.gov/cmip5/) repositories, where long cli- 
mate runs outputs are typically stored in the form of monthly averaged data. 

On a different note, the approach presented here seems promising specifically for the 
investigation of the atmosphere of exoplanets, because it allows for evaluating the error in 
the estimate of their thermodynamical properties due to the lack of high-resolution data. 

In this paper, we have worked on the post-processing of data. The analysis presented here 
should be complemented with an additional investigation of how changing the resolution of 
a model impacts the estimate of its material entropy production, in total, and process by 
process. Using the material entropy production as cost function for addressing the interplay 
between the respective role of changes in the resolution of a model and of changes in the 
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coarse graining of the post-processed data seems promising in the tantahzing quest for 
understanding what is a good model of a geophysical fluid. The results obtained here 
seem to have a much more general validity than for the specific case of the present Earth's 
climate. Therefore, we plan to extend the present analysis, by studying the combined effect 
of changes in the resolution of the model and in the effective resolution of the post-processed 
in a simpler geophysical fluid dynamical system like an Aquaplanet. 
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Appendix A: Spectral Analysis of the Impact of Coarse Graining on the Entropy 
Production for Diffusive Systems 

We wish to provide a simple outlook on how to interpret the results shown in this paper 
taking a spectral point of view. This is relevant in practical terms because many numerical 
models of geophysical fluids are actually implemented in spectral coordinates. We restrict 
ourselves to the contributions to the material entropy production coming from the presence 
of material fluxes transporting heat across temperature gradients. So, we do not consider 
here the term responsible for the dissipation of kinetic energy (which is weakly affected by 
coarse graining) nor the radiative terms contributing to the indirect formula. We consider 
the case of a much simpler simpler 3D continuum physical system, where heat transport 
obeys a generalized diffusive behavior. We make this choice not because we believe that 
the climate system is, in any real sense, diffusive, but because we wish to show that the ob- 
served dependence of the material entropy production on the coarse graining is in qualitative 
agreement with what would be obtained for a diffusive system. 

Let's assume that the contribution to the average rate of entropy production coming from 
the transport of heat due to the flow J across the temperature field T can be computed as: 

S^„, = — - / d^x / dtJ-V- = ---- / d^x / dt^—— = - — - / d^x / dt—-— 
-"* VTjy io T VTjy Jo r2 VTjy Jo T 

(Al) 
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Let's now make the simplifying diffusive-like assumption that J = —VGiT), where dG/dT > 
0, so that the flux is always opposed in verse to the temperature gradient; in the usual linear 
fiux-gradient approximation we have G{T) = kT, k > 0. We derive: 

^1 = \^^j^d'^j\t ^iVTp = ^^J^d'^^dt |W(T)p > (A2) 

where \I/(T) = J dTy^G'{T)/T'^. For sake of simplicity - but without loss of generality - 
we assume that our domain S is a parallelepiped of sides L^, Ly, and L^. Using Parseval's 
theorem, we derive that the rate of entropy production can be written as: 

p,q,r,s p,q,r,s 

where 

'^p,g,r,s = -- d^x dt^ exp[27n{p/L^x + q/Lyy + r/L,z - s/Tt)]. (A4) 

Performing a spatio-temporal coarse graining to the field ^ — )■ \& can be reframed as ap- 
plying a linear filter in Fourier space, as ^I/p,g,r,s -^ ^p,q,r,s = ^p,q,r,s'^p,q,r,s, where |$p,g,r,s| < 
IVp, g, r, s plus the usual complex conjugacy properties. Note that in our case we would 
like to able to apply the filtering to the T field and to the heating field — V ■ J and not 
to the function \1/ constructed here. Nonetheless, assuming that T is a slowly vary in does 
not change substantially across the domain (see also Eq. (14)), the conclusions are virtually 
unaltered. 

Equation A3 has the remarkable property that all of terms of the summation Sp'^^l are 
positive. We can interpret S^'^^^^^ as the entropy produced by processes occurring at the 
scales described by the indices, in this case A^ = L^/p, Xy = Ly/q, X^ = L^/r, and r = T/s. 
Therefore, indicating as S^^^ the value of the entropy production for the coarse grained 
fields, we obtain: 

^mat = 2_^ \^p ~^ K ~^ r)\^P,g,r,s\ \'^p,q,r,s\ = 2_^ \^P,g,r,s\ ^p]q,r,s — ^ maf (^5) 

p,q,r,s p,q,r,s 

In particular, we can associate a coarse graining on the scales A^^, Aj^, A^, and r (referred to 
the X-, y-, ^-directions and time, respectively) to a filter of the form ^p,q,r,s = if p > L^/Ax, 
or q > Ly/Ay, or r > L^/Az, or s > T/t and ^p,q,r,s = 1 otherwise. Slightly different ways 
of doing the coarse graining will result into different filters, which will be, nonetheless, 
asymptotically equivalent if the involved scales are the same. 
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The main conceptual point behind this result is independent of the shape of the of the 
domain of integration: the natural orthogonal expansion for atmospheric fields defined in 
an (approximately) spherical thin shell is given by spherical harmonics in for the latitudinal 
and longitudinal dependence and the usual Fourier expansion for the vertical direction. In 
the case of a thin spherical shell of thickness L^ situated at distance R from the center of 
the sphere, Eq. A3 can be rewritten as: 

C* = EE E Y.^ki+i{i+i)/R')\^>^,^^,\' (A6) 



with 



r-T 



'^n,i,m,s = ^^^ f dn [ *(2,^,0,t)exp[27rz(n/L,2-s/Tt)]*F(0,0)r. (A7) 

where Y{6, 0)™ are the usual spherical harmonics and Q refers to the solid angle. In this 
case, performing a spatio-temporal coarse graining to the field \1/ — )■ \& results in reduced 
value of the estimate of the entropy production: 

^ ' -J 

SLt = J211 E J2(^l+Kl + l)/R'mn,l,m,s\'\^n,l,n.,s\'<S^af (A8) 



It is now easy to relate the expression of |$n,z,m,sP to common averaging operation performed 
on climate data. A coarse graining on the vertical scale A^, on a temporal scale r and on a 
horizontal surface area a (or angular resolution cr/R"^) amounts to setting |$„,,i,m,sP = if 
/ > y^SnR'^/a (corresponding approximately to a triangular truncation T(K), where K is 



the integer closest to ^/8^^R^Ja), or n > L^/A^, or s > T/t and |$n,i,m,sP = 1 otherwise. 
Instead, performing zonal averages corresponds to setting, ceteris paribus, |$n,Lm,sP = if 
rra 7^ 0. 

The bottom line of the previous considerations is that adopting a coarser graining cor- 
responds to increasing the involved scales determining the spectral cutoff, if we assume a 
flux-gradient relationship which is consistent with the second law of thermodynamics (even if 
it is not the usual Fickian, linear relation), Eqs. A5 and A8 imply that as the coarse graining 
scales become larger, the estimate of the entropy production becomes smaller, because we the 
summation is performed over fewer and positive terms. This behavior is independent of the 
physical domain under consideration. Moreover, the previous considerations qualitatively 
apply - even if results are somewhat more cumbersome - if the relationship between flux and 
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gradient is more general than what previously assumed, e.g. if Ji = —diGi(T) (where the 
Einstein summation convention is not taken), under the condition that dGi/dT < Vz. 
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TABLE I. List of the main FAMOUS' parametrization routines for unresolved processes and 
their impact in term of heating rates ^^ on the various terms contributing to Smat- Codes: 
BL=Boundary Layer; AC = Atmospheric Convection: HD=Hyperdiffusion; OC=Oceanic Convec- 
tion; D=Diffusion; GW=Gravity Waves; C/E=Condensation/Evaporation: ML=Mixed Layer 
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TABLE IL Values of A 



qdir 



and A 



qind 
'-'mat 



obtained when considering the coarsest resolution 
in space (lower index: vm)^ in time (lower index: tm), or both. Only 2D horizontal averagings 
are considered here. Reference values for highest resolution data: 5*^^^ = 52.5 mWm~'^K~^ and 



^mat — 52-1 iTiWm ^K ^. All values are in units of mWm ^K ^ . 
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TABLE III. Values of A 



qdir 
'-'mat 



and A 



qind 
'-'mat 



obtained when considering the coarsest resolution 
either in horizontal direction (lower index: VM,h or vertical direction (lower index: vm,v)-, or in both 
(lower index: VM,h,v)- t is set to 1 y. Reference values for highest resolution data: S'^^t = 52.5 



mWm^^K-^ and 5^^^ -- 


= 52.1 mWm- 


"^K ^. All values are in units of mWm "^K 
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VM,v 


Mass weighted 5.3 36.5 34.5 12.8 52.1 45.6 
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FIG. 1. (a) Time-scale dependence of the total material entropy production (direct and indirect 
estimates); r ranges between 3600 s (model timestep) to 1.5 x 10^ s (50 years); (b) Differences 
between the exact material entropy production Smat and S^^^^ for 3600 s < r < 1 year (c) ; as in 



(b) but for the indirect estimate S\ 
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mat' 
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FIG. 2. (a) Material entropy production due to hydrological cycle and heat diffusion. Here L = 
50 years, (b) Material entropy production due to atmospheric small-scale temperature diffusion 
(continuous line). The material entropy production due to ocean turbulence (vertical and horizontal 
diffusion, mixed layer physics and convection is also reported, see [36] for details). We also show 
the vertical diffusion contribution (dotted-dashed line) to the total turbulent material entropy 
production. Note that the oceanic processes have a 12-hour timestep and have not been considered 
in the total entropy budget . 
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FIG. 3. Estimates of {SlZt)l (a), A 



odir 



qind 
'-^mat 



(d). The spatial 
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averaging considered here is given by longitudinal averages on horizontal surface. The x-axis 
reports r, the y-axis describes the extent A^ of the averaging in ° . 
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FIG. 4. Estimates of (5^*fj)j (a), A 
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(d). The spatial 



(b), {St^i)l (c), and A 

averaging considered here is given by areal averages along horizontal surfaces. The x-axis reports 
r, the y-axis describes the longitudinal extent IS.(j) of the coarse grained grid boxes. Details are 
given in the text. 
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FIG. 5. Estimates of {S^'at)l (a), and A 



odir 



(b) for r = 1 day, and of {S'JZ.tYv (c), and 



odir 



(d) for r = 1 year. The spatial averaging considered here is given by areal averages 
along horizontal surfaces and vertical averages along columns. The x-axis reports the numbers 
of vertical levels involved in the averaging, the y-axis describes the longitudinal extent Acf) of the 
coarse grained grid boxes. Details are given in the text. 
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FIG. 6. Estimates of (5™^j)j (a), and A 



qind 



(b) for r = 1 day, and of {S. 
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mat/v 



and 



oind 



(d) for r = 1 year. The spatial averaging considered here is given by areal averages 
along horizontal surfaces and vertical averages along columns. The x-axis reports the numbers 
of vertical levels involved in the averaging, the y-axis describes the longitudinal extent Acf) of the 
coarse grained grid boxes. Details are given in the text. 



43 



